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ABSTRACT 

X-ray observations of many clusters of galaxies reveal the presence of edges in surface brightness 
and temperature, known as "cold fronts". In relaxed clusters with cool cores, these edges have been 
interpreted as evidence for the "sloshing" of the core gas in the cluster's gravitational potential. The 
smoothness of these edges has been interpreted as evidence for the stabilizing effect of magnetic fields 
"draped" around the front surfaces. To check this hypothesis, we perform high-resolution magne- 
tohydrodynamics simulations of magnetized gas sloshing in galaxy clusters initiated by encounters 
with subclusters. We go beyond previous works on the simulation of cold fronts in a magnetized 
intracluster medium by simulating their formation in realistic, idealized mergers with high resolution 
(Ax ^ 2 kpc). Our simulations sample a parameter space of plausible initial magnetic field strengths 
and field configurations. In the simulations, we observe strong velocity shears associated with the 
cold fronts amplifying the magnetic field along the cold front surfaces, increasing the magnetic field 
strength in these layers by up to an order of magnitude, and boosting the magnetic pressure up to 
near-equipartition with thermal pressure in some cases. In these layers, the magnetic field becomes 
strong enough to stabilize the cold fronts against Kelvin- Helmholtz instabilities, resulting in sharp, 
smooth fronts as those seen in observations of real clusters. These magnetic fields also result in strong 
suppression of mixing of high and low-entropy gas in the cluster, seen in our simulations of mergers 
in the absence of a magnetic field. As a result, the heating of the core due to sloshing is very modest 
and is unable to stave off a cooling catastrophe. 

Subject headings: galaxies: clusters: general — X-rays: galaxies: clusters — methods: hydrodynamic 
simulations 



1. INTRODUCTION 

Observations with the most recent generation of X- 
ray telescopes have shown that clusters looking "re- 
laxed" upon first glance often have cool cores that are 
in a disturbed state. Many "cool-core" systems exhibit 
edges in X-ray surface brightness approximately concen- 
tric with respect to brightness peak of the cluster (e.^., 
iMazzotta et aIll2QQll : iMarkevitch eraIll2QQlL l2QQ3l ). X- 
ray spectra of these regions have revealed that in most 
cases the brighter (and therefore denser) side of the edge 
is the colder side, and hence these jumps in gas den- 
sity have been dubbed "cold fronts" (f or a detailed re- 
view see (|Markevitch fc Vikhlininll2QQ7[ )). The presence 
of one or more cold fronts is an indication of the mo- 
tion of gas in a cluster. Famous examples of cold fronts 
formed by ongoing major mergers are those in Abell 
3667 (.Vikhlinin et al.l [jpOll) and the "Bullet Cluster" 
(1E0657-56) (|Markevitch et al.l |2QQ2( ). However, cold 
fronts also occur in clusters which appear relaxed on the 
largest scales and show no obvious indications of ma- 
jor merging (e.g. Abeh 1795 and Abeh 2029). These 
fronts can arise due to "sloshing" motion of the core 
gas in the dark-matter dominated gravitational poten- 
tial. This process was studied in detaile d hyd rodynamic 
simulations in Ascasibar & Markevitch (j2QQ6[ ) (hereafter 
AM06), iTittlev fc Henrikseni (12005 ) (though with a dif- 
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ferent interpretation), ZuHone et alj (|2010f ) (hereafter 
ZMJIO), and Roediger et al. (2011). In these studies, 
the sloshing arises from encounters with small infalling 
groups or subclusters that gravitationally perturb the 
cluster core. Since the gas core is subject to ram pressure 
but the dark matter core is not, such perturbations can 
result in a separation between the gas and dark matter 
peaks, and consequently the core gas will begin to "slosh" 
back and forth in the gravitational potential well. 

In ZMJIO, who used high-resolution grid-based simu- 
lations, it was demonstrated that if the ICM is modeled 
as an unmagnetized, inviscid fluid, the initially smooth 
and sharp cold fronts that result from sloshing motions 
are quickly disrupted by Kelvin-Helmholtz instabilities 
caused by the large shearing velocities present across the 
cold fronts. Such instabilities make the fronts appear 
ragged and torn, in contrast to the fronts that are seen 
in X-ray observations, which appear smooth and sharp. 
Additionally, mixing in these simulations was very effi- 
cient, resulting in the heating of the initially cool gas 
in the core as it was mixed with hotter gas brought into 
contact by sloshing. Conversely, ZMJIO also showed that 
if the ICM is viscous, the instabilities are damped out 
and the resulting cold fronts retain their smooth, sharp 
shape. The smoothness of observed cold fronts in clus- 
ters of galaxies indicates that the ICM might possess a 
significant viscosity. An additional effect of viscosity in 
the ZMJIO simulations was to suppress mixing of gases 
of different entropies, preventing cluster cool cores from 
heating due to such mixing. 

An alternative mechanism for the stability of cold 
fronts, which has been discussed previously in the litera- 
ture, is the effect of magnetic fields. Strong observational 
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evidence points to the existence of magnetic fields per- 
meating; the cluste r volume (see ICarilli fc Taylorl I2QQ2I : 
iFerrari et al.ll2QQ"8l for recent reviews). One such line of 
evidence is sync hrotron radio emiss i on from sources suc h 
as radio halos (Feret ti et aTll2QQll : iGovoni et a T 2001), 
radio mini-halos in the cluster cool cores (iBur ns et all 
19921 : iBacchi et ah 2003; Ven turi et ar 2007: G itti et all 



20071 : IGovoni et ah. i2009i : iGiacintucci et aL. .20TTh . and 



radio relics. It is possible to estimate the magnetic field 
of the ICM under the assumption of equipartition be- 
tween the r elativistic electron population and the mag- 
netic field (|Pacholczv3 Il970f ). The magnetic field can 
also be estimated by comparing the synchrotron radio 
emission with inverse Compton hard X-ray emission from 
the same relativistic electrons, without the need to as- 
sume equipartition between the particles and the field. 
However, nonthermal X-ray emission from clusters has 
so far eluded confident detection-early reports of excess 
over the cluster thermal emission at high energies (Coma, 
Fusco-Femiano et al. 2004; A2163, Rephaeh et al. 2006) 
either allow an alternative thermal explanation or con- 
tradict the more recent upper limits from the Suzaku 
and SWIFT telescopes (e.g., Nakazawa et al. 2009; Sug- 
awara et al. 2009; Ajello et al. 2009; Wik et al. 2009). 
Nondetection of inverse Compton emission corresponds 
to lower limits on the average magnetic field strengths 
of 0.2 — 0.5 juG in radio halo regions (e.g., Ajello et al. 
2009; Wik et al. 2009) and a few /iG in relic regions 
(e.g., Finoguenov et al. 2010), broadly consistent with 
the equipartition radio estimates and Faraday rotation 
data. 

A second line of evidence for magnetic fields in clus- 
ters is that of Faraday rotation of polarized emission of 
radio sources. Rotation measure (KM) studies indicate 
that magnetic field strengths in clusters are on the order 
of a few fiGj with strengths up to tens of /iG in cluster 



cool cores (jPerlev T avlor 19 911: iTavlor fc Perle^ 
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Bonafede et al .1 120101 ). Additionally, in some clusters it 



is possible to derive RM maps, which are typically quite 
patchy, indicating that the coherence length of the clus- 
ter magnetic field is on the order of 10 kpc or less. High- 
resolution RM maps have been us ed to infer the clus- 
ter magnetic field power spectrum (Vogt fc Enfilin|[2QQ3t 



Murgia et al.l 12004; 'Vogt fc EnBhn '2005: 'Go voni et al.l 



2006; Guidetti et al. 2008; Govoni et al. 2010). These 
studies indicate that the magnetic field power spectrum 
is similar to a Kolmogorov type (Pb(^) /c~^/^), de- 
pending on the assumed value for the coherence length 
of the field fiuctuations. 

Magnetic fields oriented parallel to a shearing surface 
will sup press the growth of the Kelvin-Helmholtz in - 
stability (iLandau fc Lifshitzlll960l : IChandrasekharlfl96lf ). 
Whether or not the stability of cold fronts can be pro- 
vided by magnetic fields depends on the strength of the 
field and the orientation of the field with respect to the 
front surface. What is required is a field oriented par- 
allel to the front surface which has a magnetic pres- 
sure comparable to the kinetic energy per unit volume 
of the shearing fiow across the front. For a cold front 
in the galaxy cluster A366 7. [Vikhlinin et al.l (pOOl) and 
IVikhlinin fc Markevitch (2002) (hereafter VOl and V02) 
determined that the magnetic field strength required to 



stabilize the front is B ^ 10/iG, roughly an order of mag- 
nitude higher than the field strengths usually inferred 
from RM estimates and synchrotron diffuse radio emis- 
sion outside of the cooling core regions. Additionally, the 
indications from RM maps that the magnetic field is tan- 
gled with a small coherence length seems to argue against 
the existence of a large-scale field that could drape a cold 
front. 

However, the region surrounding a cold front is not 
a typical place in a galaxy cluster. As argued in 
V01/V02, a cold front moving through the intraclus- 
ter medium will cause the flow of the surrounding ICM 
to mo ve around it, creating a shear flow. ILvutikovl 
(|2006l ) pointed out that such flows around cold fronts 
and radio bubbles in galaxy clusters lead to "mag- 
netic draping." Provided the motion of the front is 
super- Alfvenic, a weak, tangled magnetic fleld will be 
stretched by this shear flow to produce a layer paral- 
lel to the the front surface. The magnetic fleld energy 
in this layer will be increased due to shear ampliflca- 
tion, possibly strong enough to s tabilize the front ag ainst 
Kelvin-Helmholtz instabilities (jKeshet et all |2010(). A 
number of previo us sin iulation works (e.g lAsai et al.l 
12004 120071 : iDursil 120071 : iDursi fc Pfrommed l2008f ) have 
demonstrated the stabilizing effect of magnetic flelds for 
cold fronts and AGN-blown bubbles in simplifled situa- 
tions, where cold "blobs" and hot bubbles propagated 
in simple, typically stratifled atmospheres and differ- 
ent fleld g eometri e s (e .g., uniform, tangled) were con- 
sidered. Takizawa (2008) simulated more realistic merg- 
ers with A/'-body/magnetohydrodynamics simulations, 
showing that magnetic flelds wrapped around merger 
cold fronts (such as in the Bullet Cluster) and stabilized 
them against instabilities. Here, we extend these stud- 
ies by simulating the formation of sloshing cold fronts 
produced by an encounter of a relaxed galaxy cluster 
with a small, infalling subcluster. Our aim is to deter- 
mine for realistic magnetic fleld strengths and conflgu- 
rations whether or not the magnetic flelds are amplifled 
and drape the cold fronts to a suflicient extent to sup- 
press instabilities that would otherwise grow and disrupt 
the smoothness of the fronts. 

The interplay between gas sloshing motions and mag- 
netic flelds in galaxy clusters may have other interesting 
implications. Observations of clusters with radio mini- 
halos in their cool cores revealed a spatial association 
between sloshing cold fronts and the observed radio emis- 
sion ( Mazzotta fc Giacintucci 2008). suggesting that the 
radio emission originates from the effects of the sloshing 
motions. This may come from turbulence generated by 
the sloshing motions, the ampliflcation of magnetic fleld 
strengths due to the associated shear flows (Keshet e t al.l 
Haio), or a combination of both. Another interesting 
observable effect is the large, amplifled magnetic fleld 
structures produced by sloshing may produce signatures 
in RM maps associated with the sloshing structures. Full 
3D MHD simulations of sloshing such as ours may deter- 
mine if such effects may be associated with cold fronts in 
relaxed clusters. 

This is a flrst paper in a series describing simulations 
of magnetized gas sloshing in galaxy cluster cores. The 
purposes of this work are to a) give technical details of 
the simulation setup, b) present results on the effect of 
sloshing of the gas and magnetic flelds of the cluster. 
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and c) serve as the groundwork for future studies which 
win explore additional physics and observable effects of 
sloshing. This paper is organized as follows: in Section 
2 we describe the characteristics of the simulations and 
the code. In Section 3 we describe the characteristics 
of gas sloshing in our simulations in the absence and 
presence of magnetic fields, in particular with regard to 
the effect of sloshing on the fields and the effect of the 
fields on the structure of the cold fronts. In Section 4 
we discuss the implications of these results. Finally, in 
Section 5 we summarize our results and discuss future 
developments of this work. Throughout this paper we 
assume a fiat ACDM cosmology with h = 0.7, = 0.3, 
and Qb = 0.02/^-^ 

2. SIMULATIONS 
2.1. Method 

In our simulations, we solve the ideal MHD equations. 
Written in conservation form in Gaussian units, they 
are: 



dp 

dt 



V • (/9V) = 



dt 



V • (pvv - BB) + = pg 



dt 



■ V • [^{E+p) - B(v • B)] = pg • V 



where 



E : 



SB 

'dt 



V • (vB - Bv) = 



(1) 
(2) 
(3) 
(4) 

(5) 
(6) 



where pth is the gas pressure, and e is the gas internal en- 
ergy per unit volume. For all our simulations, we assume 
an ideal gas equation of state with 7 = 5/3. 

We performed our simulations using FLASH 3, a 
parallel hydrodynamics/ A^-body astrophysical simula- 
tion code developed at the Center for Astrophysical 
Thermonuclear Fla s hes at the Univer sity of Chicago 
(jFrvxell et al.l [2QQQt iDubev et all [20091 ). FLASH uses 
adaptive mesh refinement (AMR), a technique that 
places higher resolution elements of the grid only where 
they are needed. We are interested in capturing sharp 
ICM features like shocks and cold fronts accurately, as 
well as resolving the inner cores of the cluster dark matter 
halos. It is particularly important to be able to resolve 
the grid adequately in these regions. AMR allows us to 
do so without needing to have the whole grid at the same 
resolution. 

FLASH 3 solves the equations of magnetohydrody- 
namics using a directionally u nsplit staggered mesh 
algorithm (USM; iLee Deang |2009). The USM al- 
gorithm used in FLASH 3 is based on a finite- 
volume, high-order Godunov scheme combined with 
a constrained transport method (CT), which guar- 
antees that the evolved magnetic field s atisfies the 
divergence-free condition (|Evans fc HawlevI [1988). In 
our simulations, the order of the USM algorithm cor- 
responds to the Piec e wise-P arabolic Method (PPM) of 
iColella fc Woodward! (|1984D . which is ideally suited for 



capturing shocks and contact discontinuties (such as the 
cold fronts that appear in our simulations). The USM 
solver in FLASH comes equipped with two methods for 
interpolating magnetic fields in a divergenceless manner 
from coarse to fine blocks on the AMR grid. The first 
is by simple injection of the magnetic field components 
from the coarser cells to the finer neighboring ce lls. The 
second is a higher-order method described in iBalsaral 
(|200l[ ). We find that the choice of either prescription 
does not affect our conclusions. 

The gravitational potential on the grid is set up as the 
sum of two "rigid bodies" corresponding to the contribu- 
tions to the potential from both clusters. This approach 
to the modeling the potential is used for simplicity and 
speed over solving the Poisson equation for the matter 
distribution, and is an adequate approximation for our 
purposes. It is the same approach that we used in ZMJIO 
for the resolution test; it will be described in detail in 
Section [221 

2.2. Initial Conditions 

Our initial conditions in these idealized simulations 
have been set up in a manner very similar to ZMJIO, 
with some differences which we elaborate on here. 

For the cluste r dark matter profile we have chosen a 
iHernquistI (|1990D profile: 
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27ra3 (r/a)(l + r/af 



(7) 



where Mtot and a are the mass and scale length of the DM 
halo. The Hernquist profile shares with the more com- 
monly employed Navarro, Frenk, fc White (1997, NFW) 
profile a "cuspy" inner radial dependence of the dark 
matter density, but results in simpler expressions for the 
mass, potential, and particle distribution functions. Be- 
cause we are interested in the consequences of the inter- 
action for only the central regions of the main cluster, 
the difference in the density dependence for large radii is 
unimportant. The corresponding gravitational potential 
has a particularly simple form: 



^DM(r) = 



tot 



(8) 



which for r a behaves as a point mass potential and for 
r <C a is approximately constant. For our simulations, 
the same profile shape is used for both the main cluster 
and the subcluster. 

For the gas temperature, we use a phenomenological 
formula: 

T(r) = ^V^^ (9) 
^ ^ 1 + r/a 1 + r/ac ^ ^ 

where < c < 1 is a free parameter that characterizes the 
depth of the temperature drop in the cluster center and 
ttc is the characteristic radius of that drop. This func- 
tional form can reproduce cluster temperature profiles 
of many observed relaxed galaxy clusters, which have 
a characteristic temperature drop in the center due to 
cooling. With this temperature profile, the correspond- 
ing gas density can be derived by imposing hydrostatic 
equilibrium: 



r) 



Po[l + — 



(10) 
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with exponents 



-1 -n —,P = l-n- ^7^. (11) 



a/ac 



a/ac 



We set n = 5 in order to have a constant baryon frac- 
tion at large radii, and we compute the value of po from 
the constraint Mgas/MoM = ^gas/^DM- This gas profile 
resembles those of most cool-core clusters since it contin- 
ues to increase with decreasing radius and does not have 
a flat core. From these profiles we can derive a radial 
dependence for the cluster entropy profile: 



S = kBTn 



-2/3 



1 + — 



-5/3 



1 



r/a. 



(3-2q!)/3 



(12) 

which resembles a power-law S(r) oc r^'^ over most 
of the radial ran ^e, in line with obse rvations of cool-core 
clusters (e.g.. iDonahue et alJl2QQ6f ). The initial radial 
profiles for the main cluster are given in Figure [TJ 

Our merging clusters consist of a large, "main" cluster, 
and a small infalling subcluster. They are characterized 
by the mass ratio R M1/M2, where Mi= MqR/(I^R) 
and M2 = Mo /{l-\- R) are the masses of the main cluster 
and the infalling satellite, respectively. The subcluster 
potential center starts at a distance d from the main 
cluster center, and with an initial impact parameter b. 
For all of our simulations in this parameter study, we 
choose the same subcluster mass ratio (i? = 5), distance 
d = 3 Mpc, and impact parameter b = 500 kpc. To scale 
the initial profiles for the various mass ratios of the clus- 
ters, the combinations M^/af, q, and ac^i/ai are held 
constant. For the main cluster, we chose ai = 600 kpc, 
ci = 0.17, and ac,i = 60 kpc, to resemble mass, gas den- 
sity, and temperature profiles typically observed in real 
galaxy clusters. In particular, our main cluster closely 
resembles A2029, a hot, massive, relaxe d cluster with a 
cool core t hat exhibits a cold front (e.R., Vikhlinin et al.l 
12005', '2006"; 'Ascasibar & Markevitch 2006). For ah of the 
simulations, we set up the main cluster within a cubical 
computational domain of width L = 2.4 Mpc on a side, 
with a finest cell size on our AMR grid of Ax = 2.34 kpc 
(see Appendix B for the results of a resolution test). 

The initial cluster velocities are chosen so that the total 
kinetic energy of the system is set to a fraction < K < 
1 of its potential energy, approximating the objects as 
point masses: 

E.iK-l)^ = iK-l)^^ (13) 



d 



(l + i?)2 d 



So the initial velocities in the reference frame of the cen- 
ter of mass are set to 



Vi 




(14) 



For the simulations presented in this work, we have set 
K = 1/2. 

Since both the gravitational potential of the main clus- 
ter and of the subcluster are modeled as rigid bodies, in 
order for the subcluster to fall into the main cluster and 
cause the resulting sloshing motions, it must be set up 
on a realistic trajectory. For this we choose to fix the 
center of the main cluster potential at the center of the 
domain, at rest for the duration of the simulation, and 



the center of mass of the subcluster is assumed to fall 
within this potential as a point mass. From the initial 
conditions at the beginning of the simulation the sub- 
cluster's position X and velocity v is integrated over a 
timestep At using the variable-timestep leapfrog method 
often used for integration of partic les in AT-body simula- 
tions (|Hocknev fc Eastwoodlll988[ ): 



1/2. 
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-1/2 
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with the coefficients C„ and £>„ given by 



C„ = -At" + -At"-i + - 
2 3 6 
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(At")2 



At" 



(15) 
(16) 

(17) 
(18) 

(19) 
(20) 



where n is the time index, i is the spatial index, and 
a is the particle's acceleration. By using time-centered 
velocities and stored accelerations, this method achieves 
second-order time accuracy. 

The gravitational potential at all points on the grid 
is assumed to be the sum of the respective potentials of 
the main cluster and the subcluster. The gravitational 
acceleration is then computed by finite differencing the 
potential. Since the frame we have chosen (with the main 
cluster fixed at the center) is not an inertial frame, we 
must also compute the inertial acceleration on the main 
cluster from the subcluster and add it to the acceleration 
from gravity. 

In a real cluster merger, the subcluster will be tidally 
stripped of its dark matter, growing smaller in mass, and 
(if it is bound to the main cluster) each subsequent pas- 
sage will induce a correspondingly weaker disturbance 
on the main cluster core, until the subcluster is fully 
absorbed into the main cluster. Additionally, the sub- 
cluster's orbit would be altered due to the effects of dy- 
namical friction from the surrounding dark matter of the 
main cluster. In our simplified model, the subcluster is 
not stripped of its mass, and there is no dynamical fric- 
tion, so it travels on a closed orbit which would eventu- 
ally take it on the same trajectory past the main cluster. 
This would result in a second disturbance that would be 
equal in magnitude to the original passage and be highly 
disruptive. In order to study the effects of the result- 
ing sloshing in isolation without the interference of sub- 
sequent crossings of the subcluster, the aforementioned 
trajectory is followed until the center of the subcluster 
reaches the opposite side of the box, at which point it 
is assumed to follow a const ant- velocity trajectory for 
the remainder of the simulation. As was seen in ZMJIO, 
that used a more realistic A/'-body representation of the 
clusters' DM components, the majority of the effects on 
the main cluster from the subcluster are due to the first 
core passage, so the essential features of the encounter 
are still captured. 

This gravitational potential setup is used for compu- 
tational simplicity, and is adequate for our qualitative 
study of hydrodynamic effects. We find that this pro- 
cedure for computing the trajectory of the subclusters' 
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TABLE 1 
Simulation Parameter Space 



Simulation 


Initial /3 


Initial Field Configuration 


Initial Aq (kpc) 


NoFields 


N/A 


N/A 


N/A 


BetalOO 


100 


Tangled 


43 


Beta400 


400 


Tangled 


43 


BetaieOO 


1600 


Tangled 


43 


Beta6400 


6400 


Tangled 


43 


Tight 


400 


Tangled 


15 


Loose 


400 


Tangled 


120 


Tangential 


400 


Tangential 


N/A 



orbit and the total acceleration on the gas reproduces 
well the general characteristics of the sloshing features 
seen in AM06 and ZMJIO. A forthcoming paper de- 
tailing quantitatively the similarities and differences be- 
tween the self-gravitating and the rigid potential sloshing 
setu ps supports our use of the ri gid-potential approxima- 
tion (jRoediger fc ZuHonell2Qlll ). 

Finally, it remains to set up the magnetic field of 
the cluster. For a realistic cluster magnetic field, it 
is important to satisfy a few basic conditions. The 
first is the magnitude of the field itself. Typical mea- 
surements from Faraday rotation and synchrotron radia- 
tion measurements suggest field strengths of of 1-10 /iG. 
This implies that the plasma j3 for the magnetic field 
= vIvb^ where = B'^ /^ti) is high, with typical 
values in the range 100-1000 (e.g., a field of 3 /iG at r = 
150 kpc in Abell 202 9 corresponds roughly to /3 ~ 600 
IVikhlinin et al.|[2006f ). Secondly, it is important to sat- 
isfy the constraint that V • B = 0. 

In order to have these two conditions simultane- 
ously satisfied, we use the f ollowing procedure , as i n 
iRuszkowski et all (j2003) and iRuszkowski fc OhI (|2010f ). 
A random magnetic field B(k) is set up in k-space on a 
uniform grid using independent normal random deviates 
for the real and imaginary components of the field. Thus, 
the components of the complex magnetic field in k-space 
are set up such that 

B,{y) = B^[N{ui)^iN{u2)] (21) 

By{]^) = B2[N{u:,)^iN{u^)] (22) 

B,{y) = B^[N{u^)^iN{uQ)] (23) 

where N{u) is a function of the uniformly distributed 
random variable u that returns Gaussian-distributed ran- 
dom values, and the values Bi are field amplitudes. We 
adopt a dependence of the magnetic field amplitude B{k) 
on the wavenumb e r |k| simil ar to (but not the same a s) 
lR.uszkowski et al.l (j2007[ ) and lR.uszkowski fc Qhl (|2010f ): 

B{k) cx A:-^^/^exp[-(A;/A;o)^]exp[-/ci/A;] (24) 

where /cq and ki control the exponential cutoff terms 
in the magnetic energy spectra. The cutoff at high 
wavenumber ko roughly corresponds to the coherence 
length of the magnetic field /cq = 27r/Ao (e.g., the scale 
of the observed patches in the RM maps, see ) and we 
vary this for a few of our simulations. The cutoff at low 
wavenumber ki = 27r/Ai roughly corresponds to Ai ~ 
r5oo/2 ~ 500 kpc, which is held fixed for all of our simula- 
tions. This field spectrum corresponds to a Kolmogorov 
shape for the energy spectrum (Pb(^) k~^^^) with 



cutoffs at large and small linear scales. This field is then 
Fourier transformed to yield B(x), which is rescaled to 
have an average value of y^Snp/ j3 to yield a field that has 
a pressure that scales with the gas pressure, i.e. to have 
a spatially uniform j3 for the initial field. For our simula- 
tions, we try different values of /3 to determine the effects 
of different initial field strengths, and different values of 
Ao to determine the effects of different initial field con- 
figuration. There are also predictions from simulations 
with magnetic fields and anisotropic heat conduction 
that in the cluster cool cores, the magnetic field orienta- 
tion may be preferentially tangential due to the heat-fiux 
bouyancy instabilit y (HBI, s ee, e.g.. "Parrish fc QuataertI 
2008; Bogdanovic e tall [2009; Parrish et al. 2009|). For 
this reason, it would be interesting to examine the effect 
of an initial field configuration where the field lines were 
preferentially tangential. For this purpose, we include a 
simulation that has a purely tangential initial field sim- 
ilar to that used in Bogda novic et al.. (.2009.) . described 
in spherical coordinates by 

Br = ^ (25) 

5^ = 2^0 sin l9 cos 20 (26) 

B^ = -Bo sin 20 sin 2(9 (27) 

Its magnitude was then scaled by the gas pressure. 

The entire set of simulations is summarized in Table [TJ 
Finally, for all configurations, a field consistent with V • 
B = is generated by "cleaning" the field of divergence 
terms in Fourier space via 

B(k) ^ (I-kk)B(k) (28) 

where k is the unit vector in k-space, and I is the identity 
operator. Upon transformation back to real space, the 
magnetic field satisfies V • B = 0. This field is then 
interpolated onto our AMR grid in such a way that the 
condition V • B = is maintained. 

Finally, to test the robustness of our initial model for 
the main cluster, we perform a test run with the main 
cluster kept unperturbed, with an average magnetic field 
strength near /3 = 100-200. Figure [T] shows the profiles of 
gas density, gas temperature, and gas entropy at the be- 
ginning of the simulation and at a later epoch (t = 5 Gyr 
after the beginning of the simulation), demonstrating the 
stability of these cluster quantities at all radii excepting 
the innermost couple of resolution elements (of width 
~ 2 kpc) due to force smoothing (a known numerical ef- 
fect due to the inability to resolve the gravitational force 
on scales smaller than the grid resolution). This change 
is small compared to the evolution observed in our sim- 
ulations of clusters undergoing mergers. However, the 
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Fig. 1. — Radial profiles of gas density, gas temperature, gas 
pressure, gas entropy, magnetic field strength, and plasma /3 at the 
epochs f = 0.0 Gyr and t = 5.0 Gyr for a single cluster evolved in 
isolation. 

magnetic field itself does evolve to a slightly different 
configuration, an effect we will discuss in Appendix A. 

3. RESULTS 

3.1. Sloshing of the Cool Core in an Unmagnetized Test 

Case 

First, we will briefly describe the sloshing process due 
to subcluster mergers as elucidated in Section 3 of AM06 
and Section 3 of ZMJIO. This general description is ap- 
phcable to all of our simulations, regardless of the details 
of the initial magnetic field or its presence or absence, but 
we will be referring in this description to our "control" 
simulation NoFields^ in which there is no magnetic field 
in the cluster gas. 

It is assumed that the subcluster has lost its gas due 
to ram pressure stripping from an earlier phase of the 
merger (although as the subcluster approaches the main 
cluster it begins to drag some of the cluster's ICM in a 
trailing sonic wake). This is seen in the first panel of Fig- 
ure[2l which shows slices of temperature through the clus- 
ter center in the plane of the clusters' mutual orbit. The 
core passage of the subcluster occurs at approximately 
t ^ l.S Gyr after the beginning of the simulation; each 
simulation is followed until t = 5.0 Gyr. As the subclus- 
ter approaches the main cluster's core and makes its pas- 
sage, the gas and DM peaks of the main cluster feel the 
same gravity force toward the subcluster and move to- 
gether towards it. However, the gas feels the effect of the 
ram pressure of the ambient medium. This fact becomes 
significant as the gas core is held back from the core of 
the dark matter by this pressure. After the passage of the 
core, when the direction of the gravitational force quickly 
changes, there is a rapid decline of the ram pressure. As 
a result from this change, the gas core experiences a "ram 
pressure slingshot" (jHallman & Markevitch 2004), where 
the gas that was previously held back by the ram pres- 
sure falls into the DM potential minimum and overshoots 




Fig. 2. — Slices through the gas temperature for the simulation 
with no fields at various epochs. Each panel is 500 kpc on a side. 
Major tick marks indicate 100 kpc distances. The color scale shows 
temperature in keV. The black dot marks the position of the cluster 
potential minimum. 

it. In addition to the gravitational disturbance, the wake 
trailing the subcluster transfers some of the angular mo- 
mentum from the subcluster to the core gas and also acts 
to help push the core gas out of the DM potential well. 

As the cool gas from the core climbs out of the po- 
tential minimum, it expands adiabatically. However, the 
densest, lowest-entropy gas quickly begins to sink back 
towards the potential minimum against the ram pressure 
from the surrounding ICM. Once again, as the cool gas 
falls into the potential well it overshoots it, and the pro- 
cess repeats itself on a smaller linear scale. Each time, a 
contact discontinuity ("cold front") is produced. Due to 
the angular momentum transferred from the subcluster 
by the wake, these fronts have a spiral-shaped structure. 
Throughout this process, higher-entropy gas from larger 
radii is brought into contact with the lower entropy gas 
from the core, and as these gases mix, the entropy of the 
core gas is increased (ZMJIO). The last panel of Fig- 
ure [2] shows that by the epoch t = 4.5 Gyr the initial 
temperature drop in the cluster center has disappeared. 

The resulting cold fronts in the simulation NoFields 
are qualitatively very similiar to the inviscid results from 
ZMJIO (the only difference is the way the gravitational 
potential is modeled, which in that case was computed 
from the self-gravity of the DM and gas components of 
the cluster). The fronts start off very smooth, but quickly 
develop large Kelvin- Helmholtz-driven billows, due to 
the large shear flows that are present across the front 
surfaces. The action of the instabilities on the fronts is 
to make them appear jagged, and eventually even disap- 
pear (instead of the smooth fronts seen in observations). 

3.2. The Effect of Sloshing on the Magnetic Field 

We explore our parameter space by varying two condi- 
tions of the magnetic field: the initial average strength of 
the magnetic field and the initial spatial configuration of 
the magnetic field. We will discuss the evolution of the 
magnetic field in these two sets of simulations separately. 

All of the simulations begin with a roughly uniform 
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Fig. 3. — Slices through the magnetic field strength with vectors indicating field direction in the x-y plane for the Beta6400 simulation. 
Each panel is 500 kpc on a side. Major tick marks indicate 100 kpc distances. The color scale shows magnetic field strength in iiG. 




Fig. 4. — Slices through the magnetic field strength with vectors indicating field direction in the x-y plane for the BetalGOO simulation. 
Each panel is 500 kpc on a side. Major tick marks indicate 100 kpc distances. The color scale shows magnetic field strength in /iG. 
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Fig. 7. — The alignment of cold fronts with strongly magnetized layers. Top panels: Temperature (keV, left) and magnetic field strength 
(yiiG, right) for the BetalOO simulation at the epoch t = 3.15 Gyr. Bottom panels: Temperature (keV, left) and magnetic field strength 
(/iG, right) for the Beta400 simulation at the epoch t = 3.75 Gyr. Vectors indicate magnetic field direction in the x-y plane. Each panel 
is 500 kpc on a side. 



/3 with small fluctuations. Since the measured field 
strengths in clusters of galaxies is uncertain by an or- 
der of magnitude or so, it is important to characterize 
the effect of varying the initial strength of the magnetic 
field on its subsequent evolution. The simulation set 
{BetalOO, Beta400,Betal600,Beta6400} explores the ef- 
fects of varying the initial magnetic field strength for a 
range of field strengths that covers the currently obser- 
vationally plausible interval. The initial magnetic field 
configuration for each of these simulations is the same, 
corresponding to a Kolmogorov power spectrum with 
ko = 27r/(43 kpc) and ki = 27r/(500 kpc) (see Section 

The sloshing motions drag the field around, amplify- 
ing it along shear fiows. Figures [3] through [6] illustrate 
the effect of sloshing on the magnetic field strength and 
direction over the course of time. The first panel in 
each figure shows the tangled magnetic field before the 



sloshing motions begin, with the radially averaged field 
strength highest at the center and steadily decreasing 
outward. Over the course of the next few Gyr as rep- 
resented in the following panels, the sloshing motions 
amplify the field and increase the size of the strongly 
magnetized region. These fields are also ordered on large 
scales, aligned largely along the cold front surfaces and 
other places in the domain where there are shear flows. 
In these layers the field strengths can be amplified up 
to tens of /iG. Within the envelope of the cold fronts, at 
later epochs the fields are once again tangled, though the 
field strengths have been increased by a factor of ~5-10 
over their initial values. From these figures some infer- 
ences can be made about the physical reasonableness of 
the initial conditions of the magnetic field. In particular, 
the sloshing motions in the BetalOO simulation result in 
magnetic field strengths that is over ~10 /iG out to radii 
of r ~ 150 kpc (Figure [6]), potentially in confiict with ob- 
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Fig. 9. — Slices through the plasma /3 for the simulations with 
varying initial (3 at the epoch f = 3.15 Gyr. Each panel is 500 kpc 
on a side. Major tick marks indicate 100 kpc distances. 



Fig. 11. — Slices through the plasma /3 for the simulations with 
varying initial at the epoch t = 4.5 Gyr. Each panel is 500 kpc 
on a side. Major tick marks indicate 100 kpc distances. 



ser vat ions. However, in each simulation we do not find 
any magnetic field strengths in excess of tens of /iG, im- 
plying that the fields will not increase without limit but 
will saturate. 

Figure [3 shows side-by-side examples of the tempera- 
ture and the magnetic field strength for two simulations 
at two different epochs (with magnetic field vectors over- 
laid), demonstrating the alignment of the strongly mag- 
netized layers with the cold fronts. These layers are sit- 
uated right underneath the cold front surfaces, and the 
field lines along these layers are stretched along the di- 
recti on of the fron t s, in a greement with the expectation 
froml Keshet et al.l (j2QlQ[ ). 

Figures [8] through [11] show slices of the plasma P 
through the center of the domain for the epochs t = 
2.5, 3.15, 3.75, and 4.5 Gyr after the beginning of the 



simulation, for the simulations with varying initial p. In 
each case, the shear flows amplify the magnetic field, de- 
creasing p. The degree to which (3 is decreased in these 
layers is dependent on the initial /3, however in each case 
the degree of amplification of the field energy is similar, 
on an order of magnitude or higher. 

Our default setup for the magnetic field spatial config- 
uration, as detailed in Section [221 is a tangled magnetic 
field. This field may be characterized by the slope of the 
initial power spectrum (which we take to be Kolmogro- 
rov), the cutoff at large k (/cq), and the cutoff at small 
k {ki). It is particularly instructive to examine the ef- 
fect of varying the cutoff at large k (small scales), since 
a magnetic field that is more tangled on smaller scales 
may be more difficult to amplify by shear amplification. 
The simulation set {Beta^OO, Tight, Loose, Tangential} ex- 




Fig. 12. — Slices through the plasma /3 for the simulations with 
varying initial configuration at the epoch t = 2.5 Gyr. Each panel 
is 500 kpc on a side. Major tick marks indicate 100 kpc distances. 




Fig. 13. — Slices through the plasma (3 for the simulations with 
varying initial configuration at the epoch t = 3.15 Gyr. Each panel 
is 500 kpc on a side. Major tick marks indicate 100 kpc distances. 



plores the effects of varying the initial magnetic field spa- 
tial configuration, with "Beta400" being the default con- 
figuration with ko = 27r/43 kpc~^, "Tight" correspond- 
ing to a higher ko = 27r/15 kpc~^, "Loose" corresponding 
to a lower ko = 27r/120 kpc~^, and "Tangential" corre- 
sponding to a tangentially oriented field (Equations [26]- 
[27j) . The initial average magnetic field strength for each 
of these simulations is the same, corresponding to /3 = 
400. Figures [T2l through [T5l show slices of the plasma /3 
through the center of the domain for the epochs t = 2.5, 
3.15, 3.75, and 4.5 Gyr after the beginning of the simu- 
lation. In contrast to our results for varying the initial 
/3, there is not a similarly strong dependence of the evo- 
lution of the magnetic field on its initial configuration. 
About 2 Gyr after the initial disturbance, there is little 
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Fig. 14. — Slices through the plasma /3 for the simulations with 
varying initial configuration at the epoch t = 3.75 Gyr. Each panel 
is 500 kpc on a side. Major tick marks indicate 100 kpc distances. 
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Fig. 15. — Slices through the plasma /3 for the simulations with 
varying initial configuration at the epoch t = 4.5 Gyr. Each panel 
is 500 kpc on a side. Major tick marks indicate 100 kpc distances. 



qualitative difference between the runs. The field struc- 
ture becomes similarly tangled and amplified to similar 
values. 

Within the sloshing region, the magnetic field becomes 
extremely tangled with large fiuctuations in strength. 
This is most apparent in the Tangential simulation, 
which begins with a relatively smooth tangentially ori- 
ented field but ends up with a randomly oriented field 
tangled on small scales within the sloshing region (see 
the last panels of Figures [T2lfT5]) . It appears that slosh- 
ing such as we are considering is more than enough to 
disrupt the tangential structure arising from the HBI, 
though future simulations including anisotropic heat con- 
duction will be required to confirm this (ZuHone et al. 
2011, in preparation). 

The field is most strongly amplified in the cold front 
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Fig. 16. — Initial (red lines) and final (at t = 5 Gyr, blue lines) profiles of the magnetic field strength for the simulations with varying 
initial (3. Dot-dash lines indicate the final maximum and minimum values for each radial bin. Bin widths are 10 kpc. 



layers confining tiie slosiiing region, but wiiat is the over- 
all increase in magnetic energy due to this amplification? 
Figure [16] shows the radial profile of the magnetic field 
strength, and Figure [TTl shows the radial profile of the 
plasma P of the gas, both at the beginning and the end 
of the simulations with varying initial p. The final pro- 
files have a very similar shape across the simulations. 
Out to a radius of r ~ 100 kpc, the average magnetic 
field strength is roughly constant. There is an increase in 
field strength between radii r ~ 100 — 500 kpc (the radii 
of the cold fronts at this point in the simulations), and 
after this radius the field strength begins to decline with 
radius. The main difference between the profiles is their 
magnitude compared with the initial field strength. For 
the simulations with higher initial field strengths (sim- 
ulations BetalOO and BetalOO) ^ the final average field 
strength within r ~ 100 kpc is lower than the initial av- 
erage field strength, though the field within the region of 
the cold fronts (r ~ 100 — 500 kpc) is still stronger. For 
the simulations with lower initial field strength (simula- 
tions BetalGOO and BetaG^OO)^ we find that the average 



field is amplified at most radii from its initial value. For 
the simulations with varying initial spatial field config- 
urations for P = 400, we find that the final profiles are 
very similar in magnitude and shape. 

Interestingly, the average field strengths for each sim- 
ulation with varying initial /3, and especially for simula- 
tions BetalOO and BetalOO (Figure US] and Figure Hlj), 
are all similar within a factor of ~3 (with correspond- 
ing field energies similar within an order of magnitude), 
despite the initial variation in field strength over an or- 
der of magnitude. Though sloshing increases magnetic 
field strengths at the cold fronts, that the final aver- 
age magnetic energy within the sloshing region ends up 
being the same regardless of the initial field strength 
or structure (within reasonable limits, corresponding to 
5 ~ a few jiG or P ^ several hundred). 

3.3. Cold Front Morphology in the Presence of 
Magnetic Fields 

We find that magnetic fields alter the gas dynamics of 
the sloshing cold fronts. We begin by examining the sim- 
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Fig. 17. — Initial (red lines) and final (at t = 5 Gyr, blue lines) profiles of the plasma /3 for the simulations with varying initial /3. 
Dot-dash lines indicate the final maximum and minimum values for each radial bin. Bin widths are 10 kpc. 



ulations that correspond to varying the initial plasma p. 
Figures [18] through [2T] show temperature slices through 
the cluster center for the epochs t = 2.5, 3.15, 3.75, and 
4.5 Gyr after the beginning of the simulation for differ- 
ent initial p. Two effects are noticeable. The first is 
that as the magnetic field strength is increased, the cold 
fronts appear more smooth and regular, even though the 
instabilities are not suppressed completely. This is par- 
ticularly true in the case of the fronts with the small- 
est radii (r ~ 50 kpc) at later times. These fronts are 
smooth and well-defined in the simulations BetalOO and 
BetalOO^ whereas they are more jagged and ill-defined in 
the simulations BetaG^OO and BetalGOO. The fronts at 
large radii (r ~ 150-200 kpc) in all simulations (with the 
exception of BetalOO) are strongly affected by K-H in- 
stabilities (we will discuss the reason for this in Section 
14. 2p . The second significant effect is that as the mag- 
netic field strength is increased, the cool gas in the very 
central part of the core (r ^ 50 kpc) is more resilient to 
the effects of sloshing. We will examine this effect more 



quantitatively in Section 13.41 

In contrast to the simulations where the initial mag- 
netic field strength was varied, the effect of varying the 
initial magnetic field configuration is not as significant. 
Figure [22] shows temperature slices through the cluster 
center for the epoch t = 4.5 Gyr after the beginning of 
the simulation. The smoothness of the fronts is about 
the same for each of the simulations. The temperature 
contrasts across the fronts are also very similar, and the 
temperature of the core (r ^ 50 kpc) is approximately 
the same. In conformity with the results of the simula- 
tions with varying initial /3, the cold fronts at large radii 
are most susceptible to the onset of the K-H instability. 

One other particular aspect of the temperature maps 
deserves comment. Though in many cases, the front sur- 
faces are smooth due to the large magnetic field amplifi- 
cation in their surrounding layers, the effect of the mag- 
netic pressure can become dynamically important. This 
is particularly true in the case of the BetalOO simulation, 
as in many places in this simulation, the magnetic pres- 
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Fig. 18. — Slices through the gas temperature for the simulations 
with varying initial /3 at the epoch t = 2.5 Gyr. Each panel is 500 
kpc on a side. Major tick marks indicate 100 kpc distances. The 
color scale shows temperature in keV. 




Fig. 19. — Slices through the gas temperature for the simulations 
with varying initial ^ at the epoch f = 3.15 Gyr. Each panel is 500 
kpc on a side. Major tick marks indicate 100 kpc distances. The 
color scale shows temperature in keV. 

sure is a sizeable fraction of, or even comparable to, the 
gas thermal pressure (i.e., /3 ^ 1). The result is fluc- 
tuations in density and temperature that anticorrelate 
with the fluctuations in magnetic pressure. An example 
of this is shown in Figure [23l The sum of thermal and 
magnetic pressures is relatively smooth, whereas the gas 
thermal pressure itself shows small-scale variations. The 
resulting fluctuations in density and temperature are on 
the order of ~ 10 — 30% in this simulation. 

3.4. The Effect of Magnetic Fields on Mixing and Core 
Heating Due to Sloshing 

Sloshing motions bring higher-entropy cluster gas at 
higher radii into contact with the low-entropy gas of the 




Fig. 20. — Slices through the gas temperature for the simulations 
with varying initial /5 at the epoch t = 3.75 Gyr. Each panel is 500 
kpc on a side. Major tick marks indicate 100 kpc distances. The 
color scale shows temperature in keV. 




Fig. 21. — Slices through the gas temperature for the simulations 
with varying initial /S at the epoch t = 4.5 Gyr. Each panel is 500 
kpc on a side. Major tick marks indicate 100 kpc distances. The 
color scale shows temperature in keV. 

cool core, making it possible for these gases to mix and 
the average entropy per particle of the core to increase 
(ZMJIO). We have shown in the previous paper that in 
the absence of magnetic fields, the entropy of the core 
can increase significantly from its initial state. The left 
panel of Figure [24] shows the evolution in average en- 
tropy within a radius r < 25 kpc (well within the cooling 
radius for a relaxed cluster) for each of our simulations 
with varying P and the NoFields simulation, the latter 
using the same hydrodynamical model as the invsicid 
simulations of ZMJIO. The center is taken to be the 
cluster potential minimum, and "entropy" is defined as 

2/3 

S = kBTue . Every simulation has an initial entropy 
per unit mass increase in the center associated with the 
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Fig. 22. — Slices through the gas temperature for the simulations 
with varying initial configuration at the epoch t = 4.5 Gyr. Each 
panel is 500 kpc on a side. Major tick marks indicate 100 kpc 
distances. The color scale shows temperature in keV. 



passage of the subcluster and the removal of low-entropy 
gas from the cluster core. Following this initial increase, 
the evolution of the central entropy is strongly dependent 
of the details of the magnetic field. The right panel of 
Figure [24] shows the final entropy profiles of these sim- 
ulations compared to the initial entropy profile. In the 
NoFields simulation, the cool core has been heated by 
sloshing and transformed to an isentropic core with a rel- 
atively high entropy. For the simulations with magnetic 
fields, as the initial magnetic field strength is increased, 
the final entropy profile changes less; even though slosh- 
ing brings hot and cold gases in contact, the magnetic 
field suppresses their mixing (as did isotropic viscosity 
in ZMJIO). 

In contrast, we do not find that varying the initial mag- 
netic field spatial configuration changes the effectiveness 
of mixing. Figure [25] shows the evolution in average en- 
tropy within r < 25 kpc for simulations with different 
initial field configurations (all with initial P = 400), as 
well as our control NoFields simulation. In each of the 
magnetized cases, the increase of entropy per unit mass 
of the gas in the cluster core is suppressed (in comparison 
to the NoFields simulation, but the degree of this sup- 
pression is essentially independent of the initial field con- 
figuration. This is expected, given the weak dependence 
of the final field configuration on the initial configuration 
that we noted in Section [3^ 

In this analysis so far we have ignored the effects of 
radiative cooling. In ZMJIO, a series of sloshing sim- 
ulations including radiative cooling were performed, to 
determine if the mixing of gases due to sloshing provided 
sufficient heat to the cluster core to offset a cooling catas- 
trophe. It was found in that study that sloshing was able 
to heat the core for a short period, but that this heating 
was insufficient to completely offset radiative cooling. It 
was also found that if the ICM is significantly viscous, 
mixing of hot and cold gases would be suppressed, and 
radiative cooling would be offset by a nearly insignificant 
amount of heating. Since we find that the effect of the 



magnetic field is to similarly suppress mixing and the 
heating due to mixing, we expect that a cooling catas- 
trophe would happen very quickly even in the presence 
of sloshing if the ICM is magnetized. 

To confirm this, we have performed a simula- 
tion with radiative cooling included. The cool- 
ing function is calculated usin g the MEKAL model 
(|Mewe. Kaastra. fc Liedahlll995[ ). assuming the gas has 
a uniform metallicity Z = O.SZ©. Otherwise, the simula- 
tion is identical to the Beta^OO simulation. We have fol- 
lowed the strategy of ZMJIO in beginning the simulation 
shortly after core passage (corresponding to the epoch t 
= 2.0 Gyr in the non-radiative simulations), to avoid the 
compression of the cluster core that occurs during this 
period, which would lead to significant overcooling and 
hasten a cooling catastrophe. Our aim is to determine 
the effect of sloshing on the cooling cluster core in iso- 
lation. To compare with this simulation, we have also 
ran a simulation of an isolated, magnetized galaxy clus- 
ter with radiative cooling. We follow the simulations for 
2 Gyr. 

Figure [26] shows temperature slices at the epoch t = 
3.75 Gyr for the radiative and non-radiative sloshing sim- 
ulations. The most obvious difference between the two 
simulations is the presence of the large temperature drop 
(T < 1 keV) in the center that has occurred as a result 
of the unabated cooling of the cluster core gas. The 
continued presence of low-entropy gas has also resulted 
in cold fronts with a higher temperature contrast when 
compared to the non-radiative case. Figure [27] shows 
the evolution of the average entropy with a radius of 
r = 50 kpc for the radiative and non-radiative sloshing 
simulations, compared to the radiative simulation of a 
cluster in isolation. The entropy of the core in the ra- 
diative simulation drops very quickly, and is only mod- 
estly slowed down by the heat contribution from sloshing 
when compared to the case when there is no sloshing at 
all. This shows that the presence of the magnetic field 
suppresses mixing of hot and cold gases to such a degree 
that the effect of sloshing on the cooling rate of the core 
is very modest. In line with the results from ZMJIO, we 
find that the heat delivered from sloshing alone is unable 
to stave off a cooling catastrophe, especially in the case 
of a magnetized ICM. However, thermal conduction has 
not been considered in these simulations. Since slosh- 
ing brings hot and cold gases into close contact, it may 
yet result in significant heating of the cluster core due to 
thermal conduction. An investigation of core gas sloshing 
including the effects of anisotropic thermal conduction is 
the subject of a future paper (ZuHone et al. 2011, in 
preparation). 

3.5. Mock Observations 

The observational signatures of sloshing cold fronts are 
the spiral-shaped bright edges in the X-ray emission in 
the cluster cores. To compare our simulations more di- 
rectly with these observations, we have constructed mock 
X-ray surface brightness observations from our simula- 
tion data. Each cell in the AMR grid has a photon emis- 
sion (in photons s~^ cm~^) given by 

= neUH^^iT.Z) (29) 
where rig and uh are the electron and hydrogen 
densities, respectively, and A^(T, Z) is the emissiv- 
ity which depends on temperature and metallicity. 
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Fig. 23. — An example of temperature and density fluctuations as a result of magnetic pressure from the BetalOO simulation at the 
epoch f = 3.15 Gyr. Top panels: Gas density, temperature, and magnetic pressure slices through the center of the domain. Bottom panels: 
Proflles of the density, temperature, thermal pressure, and total (magnetic + thermal) pressure along the black lines in the top panels. 
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Fig. 24. — The effect of magnetic flelds of varying initial (3 on the entropy of the gas. Left: Evolution of the average entropy within 
a radius of 25 kpc for the simulations with varying initial (3, compared to the simulation with no magnetic flelds. Right: Final entropy 
proflles (at t = 5 Gyr) for the simulations with varying initial /3, compared to the simulation with no magnetic flelds. 
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Fig. 25. — Evolution of the average entropy within a radius of 
25 kpc for the simulations with varying initial field configuration, 
compared to the simulation with no magnetic fields. 

which are assumed constant over one FLASH ceh size. 
The emissivity is calculate d using the MEKAL model 
(|Mewe. Kaastra. fc LiedahUflQQ S). under the assumption 
that the cluster is situated at redshift z = 0.06 and a con- 
stant metallicity of Z = O.SZq, an assumption adequate 
for our qualitative comparisons. Using this relation, the 
photon luminosity of each sky pixel in photons s~^ in- 
tegrated over the chosen energy range (in this case, the 
Chandra band of 0.5-7.0 keV in the observer's frame) and 
projected along the line of sight is given by 



X 



ydV « ^ Aj^ine,inH,iAVi 



(30) 



where the subscripts i refer to the quantities in each 
AMR cell. The resulting surface brightness map is uni- 
formly gridded to the resolution of the smallest cell in 
the simulation. Ax ^ 2 kpc. 

Figures [28l through [311 show the resulting X-ray surface 
brightness maps for the simulations with varying /3 for 
the epochs t = 2.5, 3.15, 3.75, and 4.5 Gyr after the be- 
ginning of the simulation. There is a trend of increasing 
smoothness of cold fronts with increasing magnetic field. 
Additionally, the simulations with higher magnetic field 
have a brighter central core, in keeping with the result 
from Section [3l4l that the cool core of the cluster is main- 
tained in these simulations. These results are similar to 
the results that were obtained in ZMJIO when viscosity 
was implemented. A more quantitative comparison with 
observations will be given in a future paper. 

4. DISCUSSION 

4.1. A Comparison of Amplified Field Strengths with 
Earlier Estimates 

Under the assumptions of ideal MHD, the magnetic 
field lines are "frozen" into the flow, a condition that 
can be expressed by the equation ( Chandr asekhar.. 196 1): 

— + V ■ (vB - Bv) = (31) 
which, combined with the continuity equation, gives 



which implies that along shear flows the magnetic field 
will be stretched and amplifled. iKeshet et al.r(|2010l ) de- 
rived an analytic estimate for the shear ampliflcation of 
the magnetic fleld along cold fronts. From Equation 11 
from Keshet et al. (2010), we have (under the assump- 
tion of compressibility): 



B(f) 

Brp' 



tdr'V - lOMiT/ 



1/2 



10 kpc 



10^ yr 



(33) 



where Mi is the Mach number of the flow inside the front, 
T4 = T/(4 keV) is the temperature inside the front, A 
is the thickness of the shear layer in kpc, and t is the 
development time of the magnetization layer. 5^ is the 
amplifled fleld strength just under and parallel to the 
front, and Br' is the initial magnetic fleld perpendicu- 
lar to the front surface. To compare with the amplifl- 
cation of p seen in our simulations, we flrst note that 
for an initially random, tangled fleld the magnetic en- 
ergy in the component perpendicular to the front surface 
is roughly one-third of the magnetic energy in the total 
fleld, 52 /Stt- (1/3)^2 

we expect that B'^/Stt 

and Bf = B^jy^ assuming B'i ^ and rearranging in 



TT. For the flnal amplifled fleld, 
TT. Setting Bi = \^Br' 



terms of /3, we flnd: 
pf 0.03 



B^l 



Pi MfT^ VlOkpcy \10^ yr 



(34) 



where f3i and Pf are the initial and flnal plasma /3. Tak- 
ing representative numbers from our simulations, we have 
Mi - 0.3 - 0.5 for the Mach number and T - 4 keV for 
the temperature just inside the cold front. In our simu- 
lations the thickness of the shear layer is A ~ 5 — 10 kpc, 
and the amplifled magnetic layers typically take a few 
X 10^ yr to develop. Under these conditions, the decrease 
in P should be Pf/Pi ~ 0.1 — 0.01, correspondng to a 
magnetic fleld energy ampliflcation B'j/Bf ~ 10 — 100, 
which is agreement with wh at we flnd in our sini ulations, 
corroborating the results of lKeshet et al.l (|2010f ). 

4.2. Smooth and Sharp Cold Fronts with Magnetic 
Fields 

Our simulations demonstrate that for magnetic fleld 
strengths compatible with those inferred from obser- 
vations, sloshing will result in shear ampliflcation and 
stretching of the magnetic fleld lines along the cold fronts. 
The stretched and amplifled flelds will suppress instabil- 
ities of the fronts and help preserve the fronts' smooth 
shapes. The efficacy of this is somewhat dependent on 
the initial strength of the magnetic fleld, as stronger ini- 
tial flelds are amplifled more quickly to strengths that 
can suppress the instabilities. 

In addition to suppressing these instabilities, the mag- 
netic flelds have another related effect which results in 
the persistence of smooth and high-contrast cold fronts. 
By suppressing the mixing of high and low-entropy gases, 
the magnetic fields ensure that the densest gas in the 
cluster remains cold c ompared to the surrounding gas. 
iGhizzardi et al.l (|201 0*) pointed out a correlation between 
clusters with sloshing cold fronts and clusters with steep 
entropy gradients in the core. This was predicted in 
AM06; if the entropy gradient in the cluster core is not 
significant, the entropy contrast is not sufficient for the 
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Fig. 26. — Slices through the gas temperature for non-radiative and radiative simulations of gas sloshing at the epoch t = 3.0 Gyr. The 
color scale shows the temperature in keV. Each panel is 500 kpc on a side. 




Fig. 27. — Evolution of the average entropy within a radius of 
50 kpc for simulations with and without radiative cooling. 



cool gas pushed out of the disturbed cluster potential 
minimum to flow back and for sloshing to begin. Simi- 
larly, sloshing will persist in generating cold fronts with 
higher contrast if the steep entropy gradient is main- 
tained during the sloshing period. Since sloshing mixes 
high and low-entropy gas, possibly eliminating this gra- 
dient (ZMJIO), mixing should be suppressed in order 
to maintain the sloshing. Our simulations indicate that 
the stronger the magnetic field is, the longer the original 
steep entropy gradient is maintained. 

In V01/V02, the lack of evidence for the growth of the 
Kelvin-Helmholtz instability on the surface of a promi- 
nent cold front in the merging galaxy cluster A3667 was 
used to argue for the existence of a magnetic field of B ^ 
10 /iG parallel to the front. This was on the basis of a 
simple stability analysis for tangential perturbations on 
a shearing surface in the presence of a magnetic field. It 
is instructive to see how the same analysis fares when ap- 



plied to our simulated cold fronts. In particular, can the 
magnitude of the actual magnetic field strength in our 
simulation be predicted? For this purpose we have cho- 
sen a few fronts in our simulations and have carried out 
the same rough calculation, made more straightforward 
by the fact that we have direct access to the relevant 
quantities (density, velocity, etc.). 

The dispersion equation for small tangential- 
discontinuity perturbations in a perfectly con- 
ducting, incompressible plasma can be written as 
(|Landau fc LifshitzlllQGOl in Gaussian units) 



47r 



47r 



(35) 



where Bh and Be are the magnetic field strengths in the 
hot and cold gases, respectively, ph and pc are the gas 
densities in the hot and cold gases, v is the shear velocity 
difference across the front, and uj and k are the pertur- 
bation frequency an d wave number. Th e discontinuity is 
stable if (VOl, VQ2. lKeshet et al.ll2QlQf ) 

5? 5? 
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(36) 



The treatment of the gas as incompressible, which sim- 
plifies the analysis, is justified by the fact that the Mach 
numbers of the gas fiows in these regions are relatively 
low (M<0.5), and the the growing modes of the K-H 
instability in the cool gas have low phase speed. 

We examine two fronts from our simulations with vary- 
ing initial /3, one that appears at the epoch t = 2.5 Gyr 
and another that appears at the epoch t = 3.15 Gyr. 
These fronts have been singled out due to their smooth 
appearance in the simulations with higher magnetic field, 
but the evidence of instability when the magnetic field 
strength is lower. The two fronts and the front cross- 
sections along which we examine the gas quantities are 
shown in Figure [321 In the simulations BetalOO and 
Beta^OO^ the chosen fronts are smooth and sharp, but 
in the simulations BetalGOO and BetaG^OO the fronts are 
visibly disturbed by instabilities. 
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TABLE 2 
Cold Front Stabillity Analysis 



Front 


Av (km/s) 


Pc/Ph (10-26 g cm-3) 


^pred (pG) 
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250 
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13.1 
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6.5/3.5 


8.0 
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Fig. 28. — Projected X-ray brightness (in the 0.5-7.0 keV band) 
for the simulations with varying (3 at the epoch t = 2.5 Gyr. The 
brightness scale is square root and the same for each panel. Each 
panel is 500 kpc on a side. 



Fig. 30. — Projected X-ray brightness (in the 0.5-7.0 keV band) 
for the simulations with varying $ at the epoch t = 3.75 Gyr. The 
brightness scale is square root and the same for each panel. Each 
panel is 500 kpc on a side. 





Fig. 29. — Projected X-ray brightness (in the 0.5-7.0 keV band) 
for the simulations with varying /3 at the epoch t = 3.15 Gyr. The 
brightness scale is square root and the same for each panel. Each 
panel is 500 kpc on a side. 



Fig. 31. — Projected X-ray brightness (in the 0.5-7.0 keV band) 
for the simulations with varying f3 at the epoch t = 4.5 Gyr. The 
brightness scale is square root and the same for each panel. Each 
panel is 500 kpc on a side. 
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For both fronts, for simplicity we have assumed that 
the densities across each front and the velocity shear are 
the same for all four simulations with different initial /3, 
which is accurate to approximately 10% and is sufficient 
for our current purpose. Table El shows the densities, 
shear velocities, and predicted minimum magnetic field 
strengths for the two fronts if they are stable, and the ac- 
tual total magnetic field strengths in the four simulations 
BetalOO,Beta400,Betal600 and BetaG^OO. 

In the BetalOO simulation, the magnetic fields at the 
front surfaces are stronger than the minimum value re- 
quired for stability, consistent with the stability of the 
fronts. In the case of the Beta^OO simulation, the fronts 
are stable, while the actual magnetic field strengths are 
very close to the values necessary to stabilize the front. 
In the BetalGOO and BetaG^OO simulations, the fronts 
are not stable, and the corresponding magnetic field 
strengths are far less than the values required for sta- 
bility. 

Strictly speaking, this analysis is only valid in the 
incompressible limit for a plane-parallel surface (see 
IChurazov fc Inogamovll2004l . for an alternative hypothe- 
sis for front stability based on the front curvature). Given 
the approximations of the above qualitative stability esti- 
mate, the predicted field estimates based on this stability 
analysis agree quite well with the actual field strengths 
that are capable of stabilizing the fronts. 

4.3. Amplification of Magnetic Fields by Sloshing: 
Implications 

Magnetic fields amplified by sloshing motions may have 
other observable effects. One such effect that has been 
suggested by previous works is that of radio mini-halos. 
Mini-halos are regions of diffuse synchrotron emission 
found in the cooling core regions of some relaxed clus- 
ters. A possible c onnection to mini-halos a nd s loshing 
was suggested by iMazzotta fc Giacintuccil (j2008l ). who 
discovered a correlation between the radio minihalo emis- 
sion and the regions bounded by the sloshing cold fronts 
in two galaxy clusters. Other clusters with evidence 
of sloshing cold fr onts also sport minihal os, including 
RXC J1504. 1-0248 (jGiacintucci et al.ll201lD and Perseus 
(ZuHone et al. 2011, in preparation). 

The existence of radio mini-halos necessitates a 
source for the relativist ic electrons. Tw o models 
have b een prop osed: the hadronic model (|Dennisonl 
19801: rV^randl Il982t iBlasi fc Colafrancescd 119991: 
Dolag fc EniBlinl 120001 : iPfrommer fc En i31in' '200l aiid 
the reacceleration model ( iCassano & BrunettfllQOll; 
ICassano et al.l 120071 : iBrunetti fc Lazariani >2011). In the 
hadronic model, relativistic electrons are produced as 
byproducts of hadronic interactions of cosmic ray pro- 
tons with thermal protons in the ICM, the resulting elec- 
trons producing synchrotron radio emission. The coin- 
cidence between the radio emission and the cold fronts 
in the hadronic model scenario would be due to the 
shear- amplified magnetic fields that are produced asare- 
sult of the sloshing, as suggested by Keshet et aTl (|2010l ) 
and seen in our simulations. Alternatively, MHD tur- 
bulence driven by the sloshing motions may accelerate 
relativistic electrons via damping of magnetosonic waves 
(|Eilekl 119791 : iBrunetti fc Laza rianI 12007 ). Numerical ex- 
periments to explore the relationship between gas slosh- 
ing in relaxed clusters and radio mini-halos is the subject 



of a forthcoming paper (ZuHone et al. 2011, in prepara- 
tion). 

Since the sloshing motions in our simulations create 
strong magnetic fields ordered on large scales, it may 
be possible to detect these fields via Faraday rotation 
measurements. A field line directed along our line of 
sight will produce a rotation measure of the polarized 
radio emission given by 

nL 

RM[radm ^] = 812 / neB\\dl (37) 

where Ue is the electron number density in cm~^, B\\ 
is the parallel magnetic field strength in /iG, and L is 
the length of the source in kpc along the line of sight. 
If the field lines in the amplified layers are partially di- 
rected along the line of sight, it should be possible to 
detect them in rotation measure observations of back- 
ground radio galaxies. In Figure [33l we give examples 
of simulated RM maps from the Beta^OO simulation at 
three different epochs. In the early stages of sloshing 
(left and center panels), long, coherent structures in the 
RM map with large values are coincident with the places 
where the cold fronts can be seen in X-rays. At later 
times, the layers are not as prominent, as the field am- 
plification is weaker and the field within the cold fronts 
is more tangled (right panel). The values of the rota- 
tion measure in these maps are comp arable to those seen 
in cool-core cl usters, such as Perseus (|Tavlor et al.ll2006[ ) 
and Hydra A (|Tavlor fc Per levlll 9931 ). The spatial coinci- 
dence of a background radio galaxy with a cold front may 
be a prime opportunity to find some observational con- 
firmation of a strong magnetic layer along a cold front. 
Alternatively, the CMB its elf may be u sed as the source 
of polarized photons ( Ohno et al.ll2003l ). possibly allow- 
ing estimates of the magnetic field strength to be made 
over a large spatial area of the cluster. Due to the very 
specific alignment of magnetic field layers with sloshing 
cold fronts, and the results of this work suggest that clus- 
ters with observed cold fronts would be good candidates 
for such a study. 

Caution should be taken, however, when interpret- 
ing these results. The accuracy of magnetic field es- 
timates from the RM method depends on the statisti- 
cal independence of the fluctuations in magnetic field 
and the thermal electron density. Anticorrelated fluctua- 
tions between the electron density and the magnetic field 
strength will result in an underestimate of the magnetic 
field in these reg:i ons, since lower- d ensity fluctuations will 
be weighed less (|Beck et al.ll2003f ). This could be a pos- 
sible concern for detecting the magnetized layers seen in 
our simulations, as Figure [23] shows that they tend to 
be associated with fluctuations of underdense gas. The 
fluctuations that we see are on the order of ~10%, so this 
should not be a strong effect. In any case, the magnetic 
fields in these layers will likely be the strongest along 
the line of sight, so they should still be easily detectable 
despite this possible systematic effect. 

5. SUMMARY 

We used high-resolution magnetohydrodynamic simu- 
lations of gas sloshing in the cluster cool cores initiated 
by the infall of subclusters. We study the effect of such 
sloshing on the magnetic field in the intracluster medium 
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Fig. 32. — Fronts chosen for stability analysis. Plots are of gas density, where the scale is the same for all panels. The fronts chosen are 
marked with a green line. Left set of panels: Front 1, for the simulations with varying /3 at the epoch t = 2.5 Gyr. Right set of panels: 
Front 2, for the epoch t = 3.15 Gyr. Each panel is 150 kpc on a side. The color scale is different from that in Figures [28l through 1311 
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33. — Simulated rotation measure maps for selected epochs in the Beta^OO simulation. Each panel is 500 kpc on a side. Color scale 



shows rotation measure in rad/m 

and the effect of this field on the cold fronts seen in X- 
ray observations of relaxed galaxy clusters. We explored 
a range of initial conditions for the magnetic field, includ- 
ing varying its initial strength and spatial configuration. 

Our results show that as a result of the shear fiows 
accompanying the sloshing motions, the magnetic field 
energy is increased significantly along the cold front sur- 
faces. The degree of amplification of the magnetic field 
strength along these surfaces is up to an order of mag- 
nitude, resulting in an overall energy of the field ^^/Stt 
that is amplified by up to an order of magnitude or two 
from its average value prior to the onset of sloshing. In 
particular, in the layers along the cold fronts, fields with 
initial strengths of 0.1-1 fiG may be amplified to tens of 
jiG. The final strength of the magnetic field is dependent 
on the initial field strengths, but weakly dependent on 
the initial spatial configuration of the field. Our results 
for the field strengths along the fronts are in agreement 
with previous analytic estimates (Keshet et al. 2010). 
The sloshing motions result in a magnetic field that is 



tangled on small scales within the cluster cool core, and 
ordered structures along the fronts, regardless of the ini- 
tial field configuration. Importantly, we find that re- 
gardless of the initial field strength or configuration, the 
final field strength averaged over the sloshing region is 
very similar, implying the field strength does not increase 
without Hmit due to amplification by sloshing, but satu- 
rates. 

If the increase of the magnetic field strength is high 
enough, our simulations show that perturbations that 
are able to grow due to the Kelvin-Helmholtz instabil- 
ity and disrupt the front surfaces are suppressed. The 
degree of suppression of these perturbations is highly de- 
pendent on the magnetic field strength of the cluster. 
Because the degree of amplification of the magnetic field 
is similar across the simulations, the simulations with the 
highest initial magnetic field strength are most effective 
at maintaining smooth front surfaces. There is very little 
dependence of the effectiveness of the field to preserve the 
fronts on the initial spatial configuration of the magnetic 
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field, because the final field does not remember much of 
its initial configuration. The cold fronts at high radii 
are most susceptible to the onset of the K-H instability, 
mainly due to the weaker magnetic fields at those radii. 

The ordered magnetic fields also suppress mixing in 
the intracluster medium. Sloshing in a cool-core cluster 
brings hot and cold gas in close contact. In the absence 
of these magnetic fields, these gases mix because of the 
development of instabilities, resulting in a net increase 
in entropy per gas particle in the core (ZMJIO). We find 
that as the average magnetic field strength in the clus- 
ter is increased, the mixing of gases is inhibited, and the 
sloshing-induced heating of the core is hindered. The 
simulations with the strongest magnetic fields result in 
very little change to the radial entropy profile of the clus- 
ter. When cooling is included, we find that the magnetic 
fields suppress mixing to such a degree that the heat con- 
tributed to the cluster core from sloshing is negligible. 

The interplay of sloshing motions and magnetic fields 
in relaxed clusters may have important consequences for 
future simulations and observations. The strong mag- 
netic fields, ordered on large scales, produced by sloshing 
could potentially be detected in rotation measure maps, 
if a radio galaxy happens to be located behind a front 
in the plane of the sky. Several examples exist of radio 
minihalos associated with X-ray cold fronts in relaxed 
clusters. Our simulations show that within the slosh- 
ing region, magnetic fields are amplified, which could be 



partially responsible for the corresp ondence between th e 
radio emission and the X-ray fronts (|Keshet et al.ll2QlQ[ ) . 
The sloshing motions are also likely to generate turbu- 
lence, and it is possible that this turbulence reaccelerates 
an existing population of relativistic electrons, which will 
be addressed in a future paper (ZuHone et al. 2011, in 
preparation). 
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APPENDIX 

A. RELAXATION OF THE MAGNETIC FIELD 

The initial magnetic fields in our simulations are not in an equilibrium configuration. On the one hand, the fields 
themselves are too dynamically weak to affect significantly the density and temperature structure of the gas, so we can 
assume that this structure will not change significantly until sloshing begins. However, if left to itself, the magnetic 
field of the cluster will rearrange itself to an equilibrium configuration, which may be a different configuration from 
which the simulation was started. In our default simulation setup, this process of relaxation is occuring while the 
subcluster is approaching the main cluster core and at the very beginning of the sloshing period. Since this period is 
not long, the relaxation is incomplete at the time of the passage of the subcluster, when the effects brought on by the 
encounter become much more important. To investigate the effects of this relaxation of the field on our results, we have 
run a second simulation with the same parameters as BetalOO^ except that the cluster has been allowed to evolve in 
isolation for 6 Gyr, when the magnetic field has approached an equilibrium configuration. Since this is the simulation 
with the strongest magnetic fields, we expect that the differences between the relaxed version and the original version 
will be most pronounced. We then begin the evolution of the simulation in the same fashion as the BetalOO simulation. 

Figure [34] shows the radial profile of the plasma j3 at the end of the relaxation of the single-cluster, compared with 
the radial profile of the plasma j3 before any evolution. The final profile runs from j3 ~ 250 in the center of the cluster 
to ~400 in the outskirts of the cluster. We might expect, based on our results from the simulations with varying 
initial /3, that the effects caused by the evolution of this initial field due to sloshing would be intermediate between the 
BetalOO and Beta^OO simulations. Figure [35] shows a comparison between temperature slices of the relaxed simulation 
and the BetalOO simulation. In the relaxed- field simulation, the Kelvin- Helmholtz instabilities are still suppressed, but 
not as strongly as they are in the unrelaxed simulation. Indeed, the t = 3.75 Gyr snapshot of the relaxed simulation 
in Figure [35] looks very similar to the one for the BetalOO simulation shown in Figure [20] Finally, Figure [36] shows a 
comparison between the evolution of the average entropy of the cool core within a radius of 25 kpc vs. time between 
the relaxed- field and BetalOO simulations. We have also included the same evolution from the BetalOO simulation for 
comparison. The increase in the entropy of the core due to mixing in the relaxed simulation is intermediate between 
the increase of the BetalOO and BetalOO simulations. 

We conclude from the results of this test that the primary differences between our BetalOO simulation and its 
"relaxed" version are that the reduced magnetic field strengths present in the latter inhibit the growth of perturbations 
on the fronts less than in the former, and that mixing and the increase of entropy in the core are similarly less inhibited. 
This is in line with the conclusions from our simulations with varying initial j3. 

B. RESOLUTION TEST 

To test the robustness of our conclusions against the effects of varying resolution, we have performed a simulation 
with the same relaxed- field, j3 = 100 setup described in the previous section, but with a finest cell size of Ax ~ 1 kpc. 
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Fig. 34. — Radial profiles of the plasma (3 for the beginning of the relaxed simulation and the Betal 00 simulation. 
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half the ceh size of the set of simulations described in this work. Increasing the resolution lowers the effective numerical 
viscosity of the simulation, allowing for smaller perturbations to be resolved and grow along the front surfaces. It also 
resolves smaller-scale fluctuations in the magnetic field, making it more difficult to maintain a coherent tangential field 
structure across the front surfaces. Both effects may affect our conclusions. 

Figures [37l through [39] show slices through the temperature (in keV) and the plasma P parameter for the epochs t 
= 2.5, 3.15, 3.75, and 4.5 Gyr after the beginning of the simulation for the relaxed- field initial setup with resolutions 
of Ax ~ 2 kpc and Ax ~ 1 kpc. The overall temperature structure, as well as the degree of field amplification, is 
very similar between the two simulations. However, two differences between the simulations are immediately apparent. 
The higher resolution-simulation allows for smaller fluctuations that "wrinkle" the front surface, and in some cases 
appear as small-scale K-H ripples (see Figure [38] for a clear example). The result is that a front that was very smooth 
in the lower-resolution simulation is not smooth in the higher-resolution case. The second difference between the two 
simulations is that the magnetic field is much more tangled in the higher-resolution case. 

Regarding the inability of the fields to suppress K-H instabilities in the higher-resolution simulation, we note that 
the places where we see the appearance of ripples in the fronts when we go to higher resolution are the places where 
even in the lower-resolution simulation we note that the magnetic layers are not particularly strong or perfectly aligned 
with the front. This inhibits their effectivness in suppressing these instabilities. In these regions the suppression (or 
lack thereof) of instabilties along the fronts appears to be a resolution effect. In all the lower-resolution simulations, 
these were the cold fronts that were most susceptible to the growth of perturbations, due to the high shear velocities 
along these fronts and the weaker magnetic field. The lower-radii fronts at both resolutions are much smoother, due 
to the stronger magnetic field strengths and the lower shear velocities across these fronts. These dense, bright fronts 
are most prominent in X-ray observations of clusters, and for them our conclusions about the stability imparted to the 
front surfaces due to magnetic fields should still be valid. However, if the ICM has a significant viscosity, its effect will 
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Fig. 35. — Slices through the gas temperature in keV for the BetalOO and relaxed simulations at the epochs t 
and t = 3.75 Gyr (bottom panels). Each panel is 500 kpc on a side. Major tick marks indicate 100 kpc distances. 



2.5 Gyr (top panels) 



be to suppress instabilities and preserve the smoothness of the cold fronts, which may make this comparison academic. 
The effect of viscosity on cold fronts in a magnetized ICM will be the subject of a future paper. 
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Fig. 37. — Slices through the temperature and plasma beta for simulations of differing resolution, at the epoch t = 2.5 Gyr. The 
simulation began with /3 = 100 and was relaxed for 6 Gyr before the encounter with the subcluster. Left: Ax = 2 kpc. Right: Ax = 1 kpc. 
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Fig. 38. — Slices through the temperature and plasma beta for simulations of differing resolution, at the epoch t = 3.15 Gyr. The 
simulation began with /3 = 100 and was relaxed for 6 Gyr before the encounter with the subcluster. Left column: Ax = 2 kpc. Right 
column: Ax = 1 kpc. 
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Fig. 39. — Slices through the temperature and plasma beta for simulations of differing resolution, at the epoch t = A.h Gyr. The 
simulation began with (3 = 100 and was relaxed for 6 Gyr before the encounter with the subcluster. Left column: Ax = 2 kpc. Right 
column: Ax = 1 kpc. 



